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We show how strong steady-state entanglement can be achieved in a three-mode optomechanical 
system (or other parametrically-coupled bosonic system) by using one of the modes as a cold reservoir 
to effectively laser-cool a delocalized Bogoliubov mode. This approach allows one to surpass the 
bound on the maximum stationary intracavity entanglement possible with a coherent two-mode 
squeezing interaction. Unlike typical dissipative entanglement schemes, the entangling dynamics 
here are most effective in a regime where the effects of the engineered reservoir cannot be described 
by a Markovian Lindblad master equation. 
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Introduction- The study of highly entangled quan- 
tum states is of interest both for fundamental reasons 
and for a myriad of applications to quantum information 
processing and quantum communication. Of particular 
fundamental interest is the possibility to entangle dis- 
tinct macroscopic objects, a task made difficult by the 
unavoidable decoherence and dissipation associated with 
such systems. Equally interesting would be the ability to 
entangle photons of very different frequencies, e.g. mi- 
crowave and optical photons. 

A promising venue for the realization of both these 
kinds of entanglement is provided by quantum optome- 
chanics, where macroscopic mechanical degrees of free- 
dom can be controlled, measured and coupled using the 
modes of an electromagnetic cavity. Recent milestones 
in this field include the ability to cavity-cool a mechan- 
ical resonator to its ground state of motion [21 |2j and 
the observation of many-photon strong coupling effects 
[3]-H. A natural setting for entanglement generation is 
a three-mode optomechanical system consisting of two 
"target" modes to be entangled, which are each coupled 
to a third "auxiliary" mode. One could either have two 
optical target modes and a mechanical auxiliary mode, 
or vice-versa; both variants have recently been achieved 
in experiment |7H9J. Several theoretical studies have de- 
scribed such schemes, using the basic idea that the aux- 
iliary mode mediates an effective (coherent) two-mode 
squeezing interaction between the two target modes |TQ] - 
[12J. However, such schemes typically yield at best a rel- 
atively small amount of intra-mode entanglement (some- 
thing which we quantify more fully below). 

In this paper, we again consider generating steady- 
state entanglement of two bosonic modes in a three-mode 
system; while we focus on an optomechanical realization, 
our ideas could also be realized using superconducting 
circuits coupled via Josephson junctions [13l [14]. Un- 
like previous works, we consider the possibility of entan- 
glement via reservoir engineering [15^: we wish to tailor 
the dissipative environment of the two target modes such 
that the dissipative dynamics relaxes the system into an 
entangled state. Such dissipative entanglement has been 




FIG. 1: Schematic of one realization of a 3- mode optomechan- 
ical system, where two cavity modes are coupled to a single 
mode of a mechanical resonator. By driving cavity 1 (2) at the 
red (blue) detuned mechanical sideband, a dissipative entan- 
glement mechanism is realized. To enhance the scheme, the 
mechanical resonator is cavity cooled and optically damped 
via a coupling to a third driven cavity mode, so that its total 
damping rate 7 is greater than the cavity damping rates Hi. 

discussed in the context of atomic systems p^H2Q] and 
has even been realized experimentally [21]. The dissipa- 
tive entanglement scheme we describe is related to op- 
tomechanical cavity-cooling schemes [^ which have 
been used successfully to cool mechanical resonators to 
the ground state. In our case, one is not cooling a simple 
mechanical mode to the ground state, but rather a hybrid 
mode delocalized over both target modes. We show that 
via a judicious choice of coupling strengths, this cooling 
can be used to yield large amounts of time-independent 
intra-cavity entanglement, far greater than that found in 
previous studies or ever possible using a coherent para- 
metric interaction. We also show that the scheme is op- 
timally effective in a regime where the engineered reser- 
voir cannot be described by a simple Markovian master 
equation; this demonstrates that the concept of dissipa- 
tive entanglement is more general than the theoretical 
setting in which it is usually presented [T9ti2T] . 

System and normal modes- While our scheme applies 
to a general bosonic three mode system, we focus here 
on an optomechanical system where two optical or mi- 
crowave cavity modes are coupled to a single mode of a 
mechanical resonator (see Fig. 1); see the supplementary 
information [24] for a discussion of the opposite case of 
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two mechanical resonators and one cavity mode. 
Hamiltonian is: 



The 
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ai is the annihilation operator for cavity i (frequency cj^, 
damping rate Hii)^ b is the annihilation operator of the 
mechanical mode (frequency cjm, damping rate 7), and 
Qi are the optomechanical coupling strengths, i^diss de- 
scribes the dissipation of each mode, as well as the driving 
of the cavity modes. To achieve an entangling interac- 
tion, cavity 1 (2) is driven at the red (blue) sideband asso- 
ciated with the mechanical resonator: Udi = cJi —uu and 
^d2 = 0J2 -\- ojm • We work in an interaction picture 
with respect to the cavity drives, and write d = di -\- di 
where di is the classical cavity amplitude. We take 
1^1,2! ^ which allows us to linearize the optomechan- 
ical interaction in the usual way (i.e. drop interaction 
terms not enhanced by the classical cavity amplitudes). 
The linearized Hamiltonian in the rotating frame is thus 
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Here Gi = gidi (we take gi^di > without loss of gener- 
ality). We further focus on the resolved-sideband regime 
co'M ^ ^1,^2, which suppresses the effects of the non- 
resonant interactions in Hcr- The remaining interaction 
Hint in Eq. (21 ) has the basic form suitable for entangling 



di and ^2- on a heuristic level, the parametric-amplifier 
interaction (G2 term) first entangles ^2 and 6, and then 
the beam-splitter interaction {Gi term) swaps the b and 
di states, thus yielding the desired entanglement. 



Note that if one made the interactions in Eq. (21) 



non-resonant (e.g. by detuning the cavity drives from 
the sideband resonances by A), one could adiabatically 
eliminate the mechanical mode, resulting in a two- mode 

squeezing interaction ^tms — ^ (did2 + h.c.^ with A ~ 

G1G2/A [25 . Such an interaction naturally leads to en- 
tanglement, but the amount is severely limited by the 
requirement of stability A < ^^1^2/2. We quantify the 
entanglement using the standard measure of the loga- 
rithmic negativity En [24 . One finds that the maximum 
stationary intracavity entanglement due to the two-mode 
squeezing coupling is (for hzi = 1^2 and zero temperature) 
Ejsf = hi(l-\-2\/n) < ln2 ~ 0.7. Many suggested schemes 
for entanglement generation in optomechanical systems 
are limited by this stability requirement. 

In contrast, the resonant case we consider allows for an 
alternative dissipative entanglement mechanism capable 
of much larger values of E^. We will focus attention on 
the regime G2 < Gi where (for ki = f<i2) our linear sys- 
tem is always stable [24l . Defining the effective two- mode 



squeeze parameter r = arctanh(G2/Gi), we introduce de- 
localized (canonical) cavity Bogoliubov mode operators: 

= di cosh r -\- dl sinh r = S (r) diS^ (r) , 

= d\smhr -\- d2Coshr = S (r) d2S^ (r) . (4) 

Here, S (r) = exp rdid2 — h.c. is a two- mode squeezing 

operator. It thus follows that the joint vacuum of /3b 
is the two-mode squeezed state |r) = 6'(r)|0,0), where 
1 0,0) is the vacuum of di^d2- The entanglement of this 
state is simply E^ — 2r. 

In terms of these new operators, = 

ujm {b>^b ^ P\Pa — P'^bPb^ optomechanical 
interactions in Eq. (21) and (|3| take the simple form 

HcK = Gp\b^ ^h.c, (5) 



Hi, 



GP\b^h.c.. 



where G = y^^^i"^^^!- The mode completely de- 
couples from the mechanics (it is a mechanically-dark 
mode [26l[27]), while in the good-cavity limit of interest 
Hcr can be neglected, implying that the mode Pa has 
a simple beam-splitter interaction with the mechanics. 
Ho-\-Hirit is trivially diagonalized, resulting in hybridized 
modes P± = (/^A ± / V2 with energies um i G. The 

existence of three distinct eigenmodes (two hybrid, one 
dark) can be useful to understand entanglement (in par- 
ticular spectral entanglement [121 [28^ ) the case where 
the mechanical mode is driven by excessive thermal noise; 
we will discuss this in a future work. 

Here, we will instead exploit the fact that Eq. ([5| has 
exactly the form used for standard cavity cooling [22 , 23j. 
Thus, if we can couple the mechanical mode 6 to a cold 
reservoir, then the beam-splitter coupling Hint can be 
used to cool $A towards its ground state, resulting in a 
stationary entangled state. To achieve this, we will take 
the mechanical mode to be coupled to a third cavity mode 
which is used to laser cool its thermal occupancy towards 
the ground state by providing a source of cold damping, 
see Fig. 1. In what follows, we include the cooling cav- 
ity coupled to the mechanical resonator in the definition 
of its effective thermal bath; hence, the damping rate 
7 includes the large contribution of the optical damp- 
ing. Amusingly, our scheme is one of the rare examples 
where the enhanced mechanical damping rate resulting 
from cavity cooling is actually highly beneficial. 

Langevin equations & cavity cooling- To describe the 
cooling potential of Hint, we next use input-output the- 
ory to derive the Heisenberg-Langevin equations for our 
linearized system. These take the standard form: 

(^-iujM ~ I) ^ ~ ^ (Gidi + G24) - Vr^in, 



dt 



d 
dt 



di 
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where (ii,in,^in describes operator- valued white noise 
driving the cavity and mechanical modes, and we have 
taken the good-cavity limit ujm ^ ^ (allowing us to drop 
terms due to Hcr)- Eqs. (p| are readily solved to find 
the steady-state occupancy and correlation of the Bogoli- 
ubov modes. In the following analytic expressions, we 
take tzi = K2 for simplicity and focus on the good-cavity 
limit (though Fig. [2] includes corrections due to Hcr)- 

Imagine first that the optomechanical interactions van- 
ished, i.e. i^int = 0, and consider the behaviour of Pa-, Pb 
(defined for a fixed r > 0). Even at zero temperature, 
Pa and j3B will have a non-zero occupancy: the Bogoli- 
ubov transformation of Eq. Q implies that vacuum noise 
driving the cavities acts as effective thermal noise for 
Pa/Pb- Writing these intrinsic (^^int = 0) occupancies as 



^ = nth,j we have 



^th,A/B = ^th,i/2 cosh^ r + (nth,2/i + l) sinh^ r. 



(7) 



Here, nth,i (^th,2) represents the temperature of the ther- 
mal bath coupled to mechanical resonator 1 (2). As one 
increases the squeeze parameter r, the effective heating of 
the j3j modes becomes exponentially large, implying the 
state of the system is far from being an ideal two-mode 
squeezed vacuum state; the state is not entangled. 

Including now the effects of i^int (and taking r = 
arctanh(G2/Gi)), the dark-mode /3b is unaffected, 
whereas the occupancy of j3A is modified to 
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where nth,M represents the temperature of the mechan- 
ical bath (which includes the cooling cavity), and the 
effective "cold damping rate" of 13a by the mechanics is 
Fopt = 4G^/7. This is the familiar equation for cavity 
cooling in the good-cavity limit, where now the mechan- 
ics plays the role of a cold reservoir. For nth,M = and 
weak coupling (Fopt <^ 7)5 the Pa mode is cooled by a 
factor i<i/{n ^ Fopt). In the strong coupling limit, the 
cooling factor saturates to a value {^/{"^ ^ n). 

Thus, while even vacuum noise tends to heat Pa, Pb 
to an exponentially-large effective temperature, the op- 
tomechanical interaction of Eq. ([5| can be used to cool 
Pa- This is insufficient to generate an ideal two- mode 
squeezed state |r) (as that requires both /3a, Pb to be 
cooled). Surprisingly, it is sufficient to generate signifi- 
cant entanglement. Using the inequality of Duan et al. 
one can show that if one cools Pa so that 



(9) 



then the two cavities must necessarily be entangled [23] . 
Further, if (Sa is in its ground state while Pb is hot 
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FIG. 2: Intra-cavity stationary entanglement (as quantified 
by log negativity ^a^) as a function of the entangling inter- 
action G2- We take nth,i = '^th,2 = 0, use a mechanical fre- 
quency of cjm = 27r X 10 MHz and damping 7 = 27r x 0.8 MHz, 
and do not make the rotating wave approximation. The 
thin red solid curve corresponds to a fixed value of Gi = 
27r X 2 MHz, and a^i = = 27r x 50 kHz. One clearly sees 
a non-monotonic dependence on G2. The solid (thick) blue 
curves and short-dashed blue curves instead correspond to 
tuning Gi to the value for each G2 that optimizes the dissipa- 
tive entangling mechanism (c.f. Eq. (10)). Solid blue (dashed 
blue) corresponds to — Hi2 — 50 kHz ( — 45 kHz, 

55 kHz) and nth,M = (nth,M = 0.3). The green dash-dot line 
represents the maximum stationary entanglement achievable 
with a coherent two-mode squeezing interaction, En — In 2. 



(i.e. given by Eq. ([7|)) then in the large r limit, the en- 
tanglement is still extremely large: ~ 2r — ln2. 

To rigorously quantify the cavity-cavity entanglement, 
we compute and discuss in what follows the log negativity 
En-, which is a function of the covariance matrix; details 
are provided in the supplementary material [24 . 

Maximizing entanglement- We now see that entan- 
glement generation is more subtle than one might ex- 
pect given the simple form of Eq. (21). In particu- 



lar, if Gi is fixed, the amount of stationary entangle- 
ment is a non-monotonic function of the entangling- 
interaction strength G2 (see Fig. [2|. The dissipative en- 
tanglement mechanism discussed here directly explains 
this behaviour, as increasing G2 has two opposing effects: 
it both increases r and the delocalization of the Bogoli- 
ubov modes (enhancing entanglement), but also increases 
the effective temperature of these modes, both by increas- 
ing the effective temperature of the cavity vacuum noise, 
and by suppressing the cavity-cooling effect (as the effec- 
tive coupling G decreases with increasing G2)- 

The maximum entanglement is achieved by carefully 
balancing the opposing tendencies described above; with- 
out this optimization, the entanglement will remain 
small. For fixed couplings G'i,G2, one can optimize the 
entanglement as a function of the mechanical damping 7. 
The maximum occurs at a non-zero dissipation strength, 
which at zero mechanical temperature and in the good- 
cavity limit is simply given by 7 = 2G. This value 
simply minimizes the occupancy of Pa, and corresponds 
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FIG. 3: Intra-cavity stationary entanglement (quantified by 
^at) versus cooperativity C2, where we use an optimized 
choice for G\ (as given by Eq. (10)), and have taken k\ — K2 



and zero temperature. The solid lines correspond to differ- 
ent choices of the damping ratio 7/^1: as indicated; increasing 
^Ik increases the amount that one can cool the delocalized 
Pa mode, and hence enhances entanglement. For large C2, 
these curves asymptote to the value in Eq. (12). The dashed 
blue line is the asymptotic expression of Eq. (11). The green 
dashed-dotted line indicates ^^a^ = ln2 (the maximum E]^ 
possible with a two-mode squeezing interaction). 



to a simple impedance matching condition (i.e. the rate 
with which the mode and mechanics exchange energy 
matches the rate at which the mechanics and its bath 
exchange energy). 

More relevant to experiment is to consider k and 
7 fixed, and optimize the entanglement over coupling 
strength. Focusing on the most interesting regime where 
the cooperativity C2 = ^2/(7^) ^ 1, and considering 
the zero temperature, good-cavity limit, we find that for 
fixed G2, the optimal G\ is given by 




^opt 

7 



C2± 
272 



1/4 



(10) 

Note that for large C2, this optimal value can easily 
correspond to a strong interaction G > a^, 7. Thus, in 
this optimal regime, the effects of our "engineered reser- 
voir" (cold mechanical resonator) on the target cavity 
modes cannot be described by a Markovian dissipator in 
a master equation; this is in stark contrast to standard 
dissipation- by-entanglement schemes. 

For the optimal value of Gi above, the entanglement 
takes simple forms in two relevant limits. If hold C2 
fixed and taking the limit 7/^^ 00, we have (at zero 
temperature): 



ET [^^2] = 2r - 21n2 



(11) 



For large C2, the entanglement is almost that of a two- 
mode squeezed vacuum (i.e. E^ = 2r). Alternatively, we 
could hold the ratio j/hz fixed and let (^2 ^ 00, we have 



T-lOpt 
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(12) 



This is the strong-interaction limit, where the 13 a mode 
hybridizes with the mechanical resonator. The maxi- 
mal cooling of Pa is consequently set by the ratio j/hc 
(c.f. Eq. ([8|). The amount of entanglement here increases 
monotonically from In 2 (the maximum possible with a 
coherent coupling) as this cooling factor is increased. 

The behaviour of the stationary entanglement versus 
coupling strength is shown in Fig. |2j where we have used 
parameters similar to those achieved in recent state-of- 
the-art experiments on microwave-circuit optomechani- 
cal systems [21 [5]. We assume that a uu = 10 MHz me- 
chanical resonator is first cavity cooled to near its ground 
state, with a final damping rate of 7 = 0.8 MHz (which is 
predominantly due to the cold optical damping used for 
the cooling). By then optimally tuning the couplings to 
the target modes Gi , G2 to optimize the dissipative en- 
tanglement mechanism (while keeping them < 2.2 MHz), 
one can obtain a relatively large En ~ 2.1. This exceed 
by an order-of-magnitude the intra-cavity entanglement 
obtained in previous studies of the same system [11 , as 
well as the the maximum of In 2 possible with a coher- 
ent two-mode squeezing interaction. If this entanglement 
was used for a teleportation experiment, the maximum 
possible fidelity would be 0.89 [30l|3T]; this reduces the 
error by a factor of three compared to what would be pos- 
sible with En = In 2. Fig. [2] also shows that large values 
of En are possible even when nth,M 7^ and hzi ^ hZ2- 

In Fig. [3] we show how the stationary entanglement 
grows to dramatically large values with C2 for an opti- 
mized choice of Gi. While the parameters needed for 
such En may be out of reach in current-generation op- 
tomechanics experiments, they may be more feasible by 
implementing a superconducting circuit realization of our 
scheme [T31IT4]. 

Conclusions- We have presented a general method for 
the dissipative generation of entanglement in a three- 
mode optomechanical system. The entanglement gener- 
ated here could be verified by measuring the covariance 
matrix of the two target cavity modes using homodyne 
techniques (see e.g. [10 ). Alternatively, one could di- 
rectly use the cavity output spectra at resonance to mea- 
sure the occupancy of the Pa mode; verifying that it vi- 
olates the Duan inequality of Eq. ([9| would also confirm 
the generation of entanglement (see [24 ). 

We thank S. Chesi and L. Tian for useful conversations. 
This work was supported by the DARPA ORCHID pro- 
gram under a grant from the AFOSR. 
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SUPPLEMENTAL INFORMATION 
System stability condition 

The linearized Heisenberg-Langevin equations for our three mode system are given in Eqs. (6) of the main text. If 
we drop the noise terms, these equations take the form v = 'M. ■ v where v = (6, di^d^) and M is a 3 x 3 matrix. For 
the system to be stable, we require that the eigenvalues of M all have a negative real part. This requirement leads 
to the well-known Routh-Hurwitz stability conditions. In the simple case tzi — the stability conditions reduce to 
the following necessary and sufficient condition: 

G'>~ (13) 

Thus, if ni = n2 and G2 < G'l, the system is always stable (regardless of the magnitude of the cavity and mechanical 
damping rates). 

In contrast, if ni 7^ ^^2, then the stability conditions are modified such that the system can be unstable at (5 = 0. 
While the general form of the conditions is somewhat unwieldy, in the interesting large cooperativity limit (i.e. Ci = 
Gl/{ni^) > 0), they reduce to: 



> C7max 



1^1 - 1^2, ■ 



(14) 



2j ^ hZi^ hZ2_ 

with C = {Gf + ^2)/(7(^i + ^^2))- Note that in the 7 1^2 hmit of interest, this condition is very sensitive to the 

sign of hzi — tZ2' On a heuristic level, damping asymmetry leads to negative damping terms in the equations of motion 
for the modes Pa^Pb- For ki < i<i2 (and G > 0), it is Pa which experiences negative damping; this is overwhelmed 
by the positive cold damping Fopt provided by the interaction with the mechanics (c.f. Eq. (8) in the main text), and 
hence there is no possibility of instability. In contrast, if hzi > hZ2^ it is the dark mode $b which experiences negative 
damping; as there is nothing to offset this, the system can become unstable. 



Definition of the logarithmic negativity 

In this paper, unless specified, we use the logarithmic negativity En to quantify the degree of entanglement; this 
quantity is a rigorous entanglement monotone, and is zero for separable states. For two-mode Gaussian states of sort 
realized by the two cavity modes di , d2 in our system, it can be calculated using the expression pL^ 

En = max[0, -ln27y"] (15) 

with 



and 



S - -detV (16) 



S = det B + det B' -2 det C. (17) 
Here V is the 4x4 covariance matrix of the two modes of interest, defined via Vjj' = ^ ^A^jA^j/ + A^j/A^j^, with 
A^j = 4' - iij) and ^= X2,P2}. Here Xi = (^di + d]^ /V2 and pi = -i (^di - d]^ j \f2 and 



di,d}- 



Sij. The 



matrix B, B' and C are 2x2 matrices related to the covariance matrix V as 

B C 

B' 



(18) 



The evaluation of entanglement needs the full information of the covariance matrix. In our main text, the relevant 
system consists of two cavity modes di and ^2, or equivalently Pa and Pb- Thus we also the correlation of Pa and 
Pb^ beside their occupancy (Eqs. (7) and (8) in the main text). The only non-zero correlator (again working in the 
good-cavity limit uu ^ ^1, ^2 where the counter-rotating terms Hcr have negligible effect) is 

PbPa) = ^ , , (nth,i + nth,2 + l)sinh2r. (19) 

/ I opt + ^ (/^ + 7) 
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Entanglement of mechanical resonators mediated by a cavity mode 



While the main text focuses on a 3-mode optomechanical system having two cavity modes coupled to a single 
mechanical mode, our entanglement-by-dissipation scheme is very general, and can be applied in principle to any set 
of three parametrically-coupled bosonic modes. In particular, it could be realized in a 3-mode optomechanical system 
where two mechanical modes are coupled to a single cavity mode. The starting Hamiltonian is now 

H = cJcava^a + Yl (^^^1^^ + {^l + ^0 ^^^) ^ ^^^^ 

i=l,2 

Here bi is the annihilation operator of the i-th mechanical resonator (frequency cj^, damping rate 7^), a is the annihi- 
lation operator of the auxiliary cavity mode (frequency cjcav, damping rate z^), and gi are the optomechanical coupling 
strengths, i^diss describes the dissipation of each mode, as well as the driving of the cavity mode. 

To achieve an entangling interaction, we will drive the cavity mode at two frequencies uj^i = ^cav — ^1 and 
^d2 = ^cav +^2 (i-e. the red(blue) sideband associated with mechanical resonator 1 (2)). We can then write the cavity 
annihilation operator as a = d{t) + de"*^^^''^, where the classical cavity amplitude d{t) = (aie*^^^ + a2e~*^2^^ ^-^uJca.^t 
is obtained from the classical equations of motion, with the sideband amplitudes ci^2 determined independently by 
the two cavity drives. We take |ai^2| ^ 1 which allows us to linearize the optomechanical interaction in the usual way 
(i.e. we drop interaction terms not enhanced by the classical cavity amplitude). Working in an interaction picture 
with respect to ^0 = 0Jib\bi + ^ ^2 + ^cavd^d, the linearized Hamiltonian in the displaced frame can be written 
H = Hint + ^CR + ^diss with 

^int = Gi (^d^bi + d\l?j + G2 [db2 + 4^^) • (21) 
Here Gi = gidi (we take di to be positive without loss of generality) and ^^cr describes counter-rotating terms 



with an explicit time-dependence. Eq. (21) takes the same form as Eq. (2) in the main text and can be used for 
entangling bi and 62 by reservoir engineering. 

As discussed in the main text, the dissipative entanglement mechanism is optimal when the mode acting as the 
engineered reservoir is both cold and has a large damping rate compared to the target modes, as this allows optimal 
cooling of the Bogoliubov mode Pa- This occurs naturally here, as the reservoir mode is a cavity mode, which given 
its higher frequency (compared to mechanical modes) will be naturally closer to the ground state, and which naturally 
has a large damping rate (i.e. its damping will be much larger than the damping rates 7 of the target mechanical 
modes, by a factor ^ 10^ to 10^ in typical optomechanical experiments). An additional benefit here is that since the 
mechanical damping rate 7 is small, the cooperativity of the system can reach a much larger value (in the strong 
coupling limit, C can reach 10^). Thus the entanglement can easily reach ~ 7 (as shown in Fig. 3), corresponding 
to teleportation fidelity of 99.9%. In contrast, the two-cavities-plus-mechanics scheme in the main text requires one 
to optically damp the mechanical damping 7 to a large value. 

While the above features are extremely attractive, the downside of this scheme comes when one considers the non- 
resonant, counter-rotating terms. For the setup discussed in the main text (two cavities, one mechanical resonator), the 
counter-rotating terms only lead to a very small heating of the coupled $a mode, much like standard optomechanical 
cavity cooling (c.f. Eq. (5) in the main text). In contrast, the situation here (two mechanical resonators, one cavity) is 
more complex. In general there are two sets of non-resonant, counter-rotating terms. The first are "diaogonal terms", 
where, e.g., the cavity drive at cjcav — oji induces Stokes scattering with mechanical resonator 1. In our rotating frame, 
these take the form: 



^CR 



= (dlgidbie-^''^''' + h,c}j + (dy2dble^'^''' + h,c}j . (22) 



We also have non-resonant "off- diagonal" terms where now, e.g., the cavity drive at cjcav — ^1 induces both Stokes 
and anti-Stokes scattering involving mechanical resonator 2. These take the form: 

^CR = at^/2 (e'^^'"^^^'^i^^ + e"'^^^+^'^'^i^^2) 

^d^29i (e'^^'^-'^'^'dbi + e^^^^+^^^^d^l) + h.c. (23) 



To understand the effects of these terms, we re-express them in terms of the Bogoluibov modes Pa^Pb using Eq. 
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(4) in the main text (with the proviso that in these equations, di should be replaced with bi). We thus have: 

H^^ = Gd/^A (coshVe-^^^^^ - sinhVe^^^^') 

+Gd/3^ sinh r cosh r (e^^^^^ - e"^^^^^) + h.c. (24) 

and 

^CR = o.lg2 sinh r + /3^ cosh + e-^^^^+^^^^d (^-^;^ sinh r + cosh r^^ 

-^a^i (e'^^'-^'^'d (/3a coshr - sinhr) + e^^^^+^^^^d (^^\ coshr - sinhr)) 

(25) 

The crucial difference here is that unlike the the case discussed in the main text (2 cavities, 1 mechanical resonator), 
the non-resonant terms when written in the Pa^Pb basis enter with factors of G1/2 coshr, Gi/2sinhr. Their effects 
will thus be exponentially larger than in the case considered in the main text, as in that case, the counter-rotating 
terms only enter with a small prefactor cx G (compare against Eq. (5) in the main text). As a result, one will need 
to be extremely far into the good-cavity limit to suppress the deleterious effects of the non-resonant terms. 

To estimate how far one needs to be in the good cavity limit in this 2-mechanical resonator setup, note that the 
terms in Hq^i^ ^cr ^^^^ contain df3A^ dPs (and their Hermitian conjugates) will cause negative damping of the Pa/Pb 
modes. If this negative damping becomes too large (i.e. larger than the intrinsic damping 7 of these modes), this can 
cause instability. Insisting that the induced negative damping is less than 7 in the large cooperativity, large-r limit 
of interest (i.e. where the entanglement can be large) leads to the requirement: 



-(-) <7 (26) 



where uj = min [2co'i , 2a;2 , Ic^i — co'2|]. Assuming the optomechanical couplings Gi has been chosen to optimize 

in the absence of counter-rotating terms and with G2 fixed (as per Eq. (10) in the main text), we can re- write this 

condition in terms of the cooperativity G2 = 02/(^^7) 



C2 



3/4 

(27) 



Thus, obtaining significant entanglement (which necessarily requires a large G2) also requires one to be extremely 
deep into the good-cavity limit, i.e. uom ^ For example, in order to achieve En ~ 5, the ratio n/Co has to be 
smaller than 10~^, a value that would be extremely challenging for current optomechanics experiments. Microwave 
cavity experiments have achieved a ratio uom/f^ ^ 50 [2 ; with such numbers, an En ^ 1.5 could be possible. This is 
already more than a factor of two greater than the maximal entanglement En = In 2 possible with a direct two-mode 
squeezing interaction. 

The bound of the Bogoliubov mode occupancy by the generalized Duan inequality 

Using the defintion of Bogoliubov mode Eq. (4) in the main text, one gets 

l + 2nA 



Vii cosh^ r + V22 sinh^ r + (G12 + G12) sinh r cosh r, (28) 



with Vii = V22 = ^|<^25^2|^ and G12 = ^|<^i,<^2|^ where {• • • } denotes anti-commutator. 

According to generalized Duan inequality [3 , any separable (non-entangled state) will necessarily satisfy 




D = l^axi + -f 2J y + [api - -p2j y > ^ + «^ (29) 
where Xi = (^di + d]^ j \f2 and pi = —i (^di — d]^ j \f2 and a is an arbitrary nonzero real number. 
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If we set = cothr, we can easily verify that 

1 + 2hA = ^ sinh 2r. (30) 

For this choice of a, using the Duan bound on D thus yields that for any separable state, the occupancy of the Pa 
mode must satisfy: 

fiA > sinh^ r. (31) 

Note from Eq. (7) in the main text that at zero temperature at zero optomechanical coupling, ua is exactly sinh^ r, 
indicating no violation of the Duan bound (consistent with no entanglement in this case). However, any amount of 
cooling (as per Eq. (8)) in the main text will then lower the value of n^, leading to a violation of the Duan bound 
and hence entanglement of the two target cavity modes. 



Measuring entanglement via the $a occupancy 

As shown in the previous section, if one can verify that the occupancy of the /3a mode has been reduced below the 
value sinh^r (where r = arctanhG2/Gi), then one has violated the Duan inequality, indicating that cavities 1 and 2 
are necessarily entangled. We now show that this mode occupancy can be directly obtained from the output spectra 
of the two cavities. Consider the simplest case where all dissipative baths are at zero temperature, and where we are 
deep in the good-cavity limit; for concreteness, also focus on the output spectrum of cavity 1. The solution to the 
Heisenberg-Langevin equation of motion for di may be written: 

dl[(^] = XlH (^-iGiXM[(^] (^-iG/^^H - ^/^bin[(^]^ - AA^^l,in) (32) 

where the susceptibilities are: 

XiH = ^7 7^ (33) 

XmM = \^ (34) 

-z((j -(Jm) +7/2 

Further, the output field from cavity 1 (ii^out is given by the standard input-output relation 

^^l,out = ^^l,in + \^di, (35) 

and its power spectral density defined as 

5iH = y"dte'-*(dX,,,(i)di,out(0)). (36) 

Note that the above equations (like the equations in the main text) are in the rotating frame set by the cavity drive 
frequencies; hence, cj = cjm corresponds to the cavity resonance. 

We see immediately that the only contribution to Si[uj] will be from the non-zero occupancy of (3a- We focus for 
simplicity on the regime where the effective coupling G is less than 7, a regime where large entanglement is possible 
but the /3a mode does not hybridize with the mechanical resonator; the entanglement in this regime is described by 
Eq. (11) of the main text. In this regime, the total damping rate of the /3a mode will he Fopt, where the effective 
cold damping rate due to the mechanics is Fopt = 4G^/7 (c.f. Eq. (8) in the main text). A straightforward calculation 
then shows that 

5iM = 16Ci^^^(/3i/3^) (37) 

where the co-operativities Gi and G are defined as Gi = G\/{n^)^ C = G^/(/^7). We thus see that the output 
spectrum of cavity 1 at resonance (i.e. the number of photons leaving the cavity at its resonance frequency) gives a 
direct measure of the occupancy of the coupled Bogoluibov mode ^a- 
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